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Abstract. We present results from a simulation of SU(2) lattice gauge theory with Nf = 4 flavors of Wilson 
fermion and non-zero quark chemical potential /i, using the same 12^ x 24 lattice, bare gauge coupling, 
and pion mass in cut-off units as a previous study [9] with Nf — 2. The string tension for Nf = 4 is found 
to be considerably smaller implying smoother gauge field configurations. Thermodynamic observables and 
order parameters for superfluidity and color deconfinement are studied, and comparisons drawn between 
the two theories. Results for quark density and pressure as functions of n are qualitatively similar for 
Nf = 2 and Nf = 4; in both cases there is evidence for a phase in which baryonic matter is simultaneously 
degenerate and confined. Results for the stress-energy tensor, however, suggest that while Nf = 2 has a 
regime where dilute matter is non-relativistic and weakly-interacting, Nf = 4 matter is relativistic and 
strongly-interacting for all values of fj. above onset. 



1 Introduction 

A first principles QCD-based understanding of the be- 
haviour of dense strongly interacting matter would be of 
tremendous value, in part because of the theoretical inter- 
est in new phases of matter exhibiting properties such as 
color superconductivity and crystallisation, and also be- 
cause a quantitatively accurate equation of state s{nq), 
p{nq) (where e is energy density, p pressure and quark 
density) is a necessary input for the relativistic equations 
of stellar structure determining the properties of compact 
stellar objects. Lattice simulations, in principle the most 
reliable systematic non-perturbative approach to QCD, 
are unfortunately unavailable for dense matter, since in 
the presence of a quark chemical potential /d biasing 
the system to have more quarks than anti-quarks, the Eu- 
clidean action is complex-valued making Monte Carlo im- 
portance sampling inoperable. Our approach to the prob- 
lem is to study QC2D, ie. Yang-Mills plus fundamental 
quarks with gauge group SU(2), which can be shown to 
have a real positive functional integral measure for an even 
number of quark flavors Nf. It is possible using entirely 
orthodox lattice QCD simulation methods to generate rep- 
resentative gluon field configurations and study thermo- 
dynamic properties of this model for arbitrary fi - QC2D 
is thus the simplest theory in which dense matter with 
long-ranged inter-quark interactions can be systematically 
studied. 

Previous work has shown that as fi is raised baryonic 
matter is induced to form in the ground state at an onset 



chemical potential /io — 



since the lightest baryon 



in the physical spectrum of QC2D is degenerate with the 



pion [T1I5JI3]. At the same point the U(1)b baryon-number 
symmetry of the model is spontaneously broken by a su- 
perfluid diquark condensate ■ Studies of the hadron 

spectrum confirm that in this regime there is a Gold- 
stone excitation formed from a superposition of qq and qq 
bound states [5,2,6 . There is also evidence for a change 
in the level ordering of tt- and p-meson states as density in- 
creases [71[S]. Subsequent simulations probing higher den- 
sities have found evidence for degenerate quark matter and 
color deconfinement [8j|9]. There have also been interest- 
ing studies of the glueball spectrum 10^ and the presence 
of gauge field fluctuations with non-trivial topology [TTl 
112] in a dense medium. 

It is, however, difficult to assess this body of work sys- 
tematically, because the results have been obtained with 
differing lattice fermion formulations (staggered, Wilson), 
differing numbers of quark flavors {Nf = 2, 4, 8), and even 
differing matter representations (fundamental, adjoint). 
Moreover, away from the continuum limit (lattice spacing 
a 0) staggered lattice fermions have global invariances 
and symmetry breaking patterns distinct from those of the 
continuum theory; this has a profound effect in theories 
such as QC2D having real or pseudoreal matter represen- 
tations, even influencing the nature of the Sign Problem 
(see eg. [T]). 

In this work we attempt to remedy the situation by 
presenting results for Nf = A obtained with the same lat- 
tice action used in previous studies with Nf = 2 [5112]. 
It is intrinsically interesting to examine the influence of 
varying the matter content on the theory's ground state, 
especially since Nf /Nc > 1 in this case, although it should 
be immediately noted that the first two coefficients of the 
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renormalisation group beta-function remain negative, im- 
plying the model remains asymptotically-free and confin- 
ing for /i, T — Q- We will see as fj, is varied that some 
features of the Nf = 4 theory are qualitatively similar 
to those found for Nf = 2, while others, particularly the 
behaviour of the stress-energy tensor T^^, appear qualita- 
tively very different. Another important difference is that 
for the same bare lattice gauge coupling /? and pion mass 
mTrCt (ie. measured in cut-off units) , the string tension aa^ 
is much smaller for Nf = 4. Strictly, since the two lattice 
models yield distinct continuum theories, this cannot be 
interpreted in terms of a finer lattice; nonetheless the re- 
sulting gauge field configurations are markedly smoother 
for Nf = 4 than for Nf = 2, implying that the influence 
of at least some lattice artifacts will be diminished. In a 
follow-up paper [12j we will demonstrate that this makes 
the identification of topological excitations much easier to 
control, allowing their behaviour in a dense medium to be 
probed. 

In Sec. [5] we will detail the lattice action used, and 
give results for the string tension measured at = 0. The 
main results for quark density nq(fi), pressure energy 
density e(/i) and trace of the stress energy tensor r^^(/i) 
are presented in Sec. [31 as well as details of both superfluid 
and Polyakov line order parameters yielding information 
on the physical nature of the matter which forms. We 
will find that the main claim of namely that there 
is a range of /i in which baryonic matter in QC2D is at 
once degenerate (ie. having a well-defined Fermi-surface) 
and confined, is supported by the results of the current 
study. Both similarities and differences between Nf = 4 
and Nf = 2 are discussed in Sec. [H 



2 Simulation Details 



The lattice action used, as in previous studies in this se- 
ries [HllS] , employs Wilson fermions interacting with gauge 
fields governed by the Wilson gauge action. In units where 
the lattice spacing a = 1, 
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where Ui,\ is the oriented product of 4 SU(2)-valued link 
fields U^x around the sides of an elementary plaquette in 
the u-X plane, and fi is the quark chemical potential. The 
fields and are SA^'c-component spinors located on the 
lattice sites; if we write S'* = ((/?, 0)*, then 



fM[U;fi] Atj75 
{~Kj^5 M[U;-^l] 



(2) 



^ Note that Nf — 4 lies below the perturbative prediction 
^ ^ for the existence of a conformal fixed point. 



with M[U; 11] the standard Wilson fermion matrix given 

by 



Mxy[U ; ^ Sxy - 



y^y,x~v 



(3) 



With the identifications ipi = ip, 'tpi = ip, tp2 — (0Cr2)*'', 
'4'2 — {CT2(j)Y^ , where C'^fj^C'^ = —7* and T2 acts on 
color, the action (HI is readily seen to be equivalent to two 
copies of an action describing a quark isodoublet ("01, '02) 
with the usual coupling to gauge fields, with in addition a 
gauge invariant scalar isoscalar diquark coupling (or Ma- 
jorana mass term) of the form 



(4) 



A diquark source j mitigates long-wavelength fluctu- 
ations and hence critical slowing down in any superfluid 
phase characterised by (<^750) 7^ 0. It explicitly breaks the 
global U(1)b symmetry (p ^ e'"<p, >-)■ e~'"0 of ([T]0. Al- 
ternative choices of diquark operator consistent with the 
Pauli Exclusion Principle are possible, but in such cases 
numerical simulations show no firm evidence for symme- 
try breaking leading to superfluidity [131 . In any case, the 
"physical" limit j — >■ is potentially as technically and 
computationally challenging as the chiral limit in vacuum 
QCD. 

It is possible to show that detA^ is real and positive 
and hence that the model can be simulated using an ortho- 
dox hybrid Monte Carlo algorithm [5]. In this first study 
with Nf = 4, our strategy is to compare with results ob- 
tained for Nf = 2, using the same 12'^ x 24 lattice and 
bare gauge coupling /? = 1.9 [S]. Since in principle the two 
theories have distinct physical properties, there is some ar- 
bitrariness in identifying an appropriate matching condi- 
tion. In QC2D the onset transition at which quark density 
first rises from zero signalling the transition from vacuum 
to baryonic matter at T = is expected at fio = jWn- [H]- 
Accordingly we have chosen to match the pion mass mea- 
sured in lattice units, ensuring that the onset happens at 
a similar value of /za in the two simulations. We varied k 
whilst keeping fixed (3 = 1.9, n = j = 0, measured the 
pion propagator using a local source and sink, and found 
that K = 0.158 yielded m,ra = 0.677(14) in good agree- 
ment with the value m^a = 0.68(1) used in the Nf = 2 
study ^9]. 

With no further freedom, in order to set the physical 
scale we next measured Wilson loops with R x T with 
all possible spatial separations R. Fitting W{R, T) to the 
form Ae~^^^^"^ for aTmin > 5 (see Fig. [T]), we use the 
Cornell potential 



V{R) = C+-+aR 
R 



(5) 



In the limit j — /i — the full global symmetry group of 
the action ^ is Sp(8). 
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Fig. 1. The static quark potential from fits to Wilson loops 
with aT > 5, together with a fit to the Cornall potential ((S)). 

to extract the string tension a. We find y/cra — 0.1096(64), 
yielding a lattice spacing a = 0.052(3)fni assuming a phys- 
ical string tension (440 MeV)^ about a factor of three 
smaller than the value a — 0T86(8)fm found for Nf — 
2 [5]. Some insight may be gained by integration of the 
renormalisation group beta-function /3(g) = —adg/da: 

with the two-loop coefficients given by 

(7) 

Fig. m plots the solution ^ for A^c = 2, 3 and Nf = 2, 4. 
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Fig. 2. Plot of the RG equation solution 0) for Nc — 2,3 and 
Nf ^2,4:. ' 

Whilst it is difficult to directly match the bare lattice cou- 
pling ^2Ncl fi to the Yang-Mills coupling g plotted along 
the horizontal axis, it can be seen that with two colors 
at the "boosted" value g « 2.9, the effect of doubling the 
number of flavors has the same dramatic effect of reducing 



a by a factor of three; in fact the constant C in ([5]) can be 
chosen so that the vertical scale of Fig. [2] matches our es- 
timates for a made using the string tension a. For N^ — 3 
with the same coupling, by contrast, changing N ^ from 2 
to 4 reduces a by just 35%. The small lattice spacing at 
N f = 4 is therefore not so surprising, and underlines that 
changing the flavor content in QC2D is not as innocent as 
might naively be thought. It also challenges our identifica- 
tion of a physical scale using the string tension (particu- 
larly since we don't have the Particle Data Group to help 
in QC2D). Taking = 440MeV yields ~ 2.6GeV 
and the temperature of the lattice T ~ 160MeV, which 
on the face of it makes comparison with the results of [9] 
difficult. 

Without data at different values of a enabling a contin- 
uum extrapolation, the best we can achieve in this inital 
study of flavor-rich dense matter is therefore a compari- 
son of Nf — 2,4 in pion-mass units, since the two simu- 
lations have matched values of LsTrij^, LtmriT^, and in the 
T — 0, j ^ limits should manifest the onset of baryonic 
matter at the same value of plo/ttIt^ (recall also that from 
both theoretical [14] and simulational [1] standpoints, the 
best way to present data obtained for ^ ^ jio with differ- 
ent quark masses uses units of ^/niT,). As a result, the 
figures plotted in Sec. [3] all start to exhibit non-trivial be- 
haviour at the same point along the /i-axis. 

We should bear in mind that at Nf = 4 the gauge 
field configurations are likely to be significantly smoother 
and hence lattice artifacts smaller. A potential worry is 
that the Nf = 4 simulations are much more susceptible 
to finite- volume artifacts and thermal effects as a result of 
the much smaller values of Lg^a and Lt\Ja. 



3 Numerical Results at /i 7^ 

The action ([Ij was used to generate between 300 - 500 
HMC trajectories of mean length 0.5 for chemical po- 
tentials /ia g [0.25,1.2], with diquark source ja — 0.04 
throughout. The molecular dynamics timestep needed to 
maintain 75% acceptance ranged from dt = 0.004 at fia = 
0.25 to dt = 0.002 at fia = 1.0. Fig. [3] shows results for 
the primary thermodynamic observable, the quark den- 
sity nq{^) = —d\nZ/dfi, for both Nf = 4 and Nf — 
2 [9]. What is actually plotted is the dimensionless ratio 
riq/n^gg, where n}gg is the quark density evaluated on the 
same 12^^ x 24 lattice with free massless quarks, ie. with 
13 = 00, K — ^. The purpose of this is to compensate for 
both finite- volume effects and lattice artifacts [8j. Note 
that for T = in the thermodynamic and continuum lim- 
its, 

^sS\^) = (8) 

characteristic of a degenerate system with Ep — kp = fi. 

The most striking feature of Fig. [3] is the qualitative 
similarity of the Nf — 2 and Nf = 4 data when ex- 
pressed in cut-off units. In both cases nq/n}g^ rises sharply 
from near the theoretical onset at /io ~ ^tItt = 0.34a~^ 



4 



Simon Hands et al.: Lattice Study of Dense Matter witli Two Colors and Four Flavors 



2.5 r 



2 - 



1.5 ■ 




We note that the Nf = 4 curve is smoother, and that 
once again the ratio p/psB is close to unity in the range 
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Fig. 3. n^jn^sB andp/p^i)^ vs. ^/m^ for Nf = 2,4. 



to a maximum at /.t/m^ « 0.6 before falling. This non- 
monotonic behaviour is predicted in chiral perturbation 
theory (xPT) [T4ll8]. which treats the matter forming at 
onset as a Bose- Einstein condensate of deeply bound scalar 
isoscalar diquark pairs. In fact beyond the maximum, xPT 
predicts the ratio to be monotonically falling, proportional 
to as /i — oo. By contrast, at a value /i = /iQ « 
0.75rn.n- the measured ratio levels off; moreover, for Nf = 4 
the ratio is close to unity over the range 0.75 ^ ^J,/mTr ^ 1-0. 
Finally, for /i > no ~ l.OmTr the ratio starts to rise again. 

In [3] the behaviour for /i e {fiQ^^jj) wa,s interpreted 
as a regime where the fermionic nature of the quarks has 
become manifest, resulting in the formation of a Fermi 
surface only weakly perturbed by diquark condensation 
(see below). In other words, the threshold value /ig corre- 
sponds to a BEC/BCS crossover. The new data at Nf = 
4 strengthens this interpretation, since the ratio is now 
so close to one; indeed nq/n^g^ seems to approach unity 
from above as a — s- 0. Since, however, care must be ex- 
ercised when extrapolating to the continuum limit using 
data from distinct theories, and also because the necessary 
J — extrapolation is still to be done, we cannot say at 
this stage whether any physical significance should be at- 
tached to the differing heights of the peaks at /i/m^ — 0.6. 
The flatness of the Nf = 4 data in this regime does how- 
ever increase our confidence in the rather ad hoc procedure 
for mitigating lattice artifacts by taking the ratio of inter- 
acting to free theories, which can only strictly be justified 
if a continuum limit is taken. 

Fig. [3] also plots the pressure p for both theories, ob- 
tained by integrating the Maxwell relation Uq = dp/dp.. 
For the lattice data we have implemented this via the for- 
mula [5] 
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"-SB 



{Pi')dpi' 



(9) 



with the integrals estimated by an extended trapezoidal 
rule. Since the data for p derive from those for Uq, most 
of the qualitative comments made above apply here also. 
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Fig. 4. {qq}/)!^ and rescaled Polyakov line vs. pb/m-n for Nf 
2,4. 



In order to characterise the properties of the different 
regions better, in Fig.|4]we plot two further important ob- 
servables, the superfluid order parameter {qq) = ^{ip^Q(j)— 

and the Polyakov line (L) = ]^tr(nfii C^O;x,t)- 
Once again, the results for Nf — 4 are qualitatively very 
similar to those found for Nf — 2. Note that (qq) is strictly 
only an order parameter in the limit j — 0; we choose to 
plot the ratio {qq) / p? since for a degenerate system in 
which superfluidity arises as a result of BCS condensa- 
tion the order parameter is expected to scale as the area 
of the Fermi surface via {qq) oc AbcsP^, where Abcs is 
the superfluid gap. Fig. 2] confirms that this is indeed the 
case for p, S {fig, fio)', the data for Nf = 4 being flatter 
than those for Nf — 2, although both manifest a shallow 
minimum at fi/rriTr ~ 0.95. This suggests the inequality 



^BCS 



^BCS 



(10) 



Whilst this may signal a physical difference between the 
two theories, the influence of lattice artifacts, or indeed the 
opening up of other more favoured condensation channels 
in the flavor-rich case, cannot yet be eliminated. 

Fig. m also shows a sharp transition between (i) « 
and {L) > at fiD ~ I.Otott- This transition approxi- 
mately coincides with the upturn in nq/n}^g observed in 
Fig. [31 Since (i) ^ exp{—fq/T), where fq is the free en- 
ergy of a static isolated color source, we interpret pu as 
the chemical potential at which color deconfinement takes 
place. It is striking that fj,^^~'^ /ttIt^ and //^^"^/tott appear 
to be identical; we attribute the difference in magnitude of 
(L) for > /z_D to the smoother gauge fields for Nf = 4, 
resulting in a smaller downwards renormalisation of the 
Polyakov line [iQ . 

We therefore recover for QC2D with Nf = 4 the same 
intriguing result found for Nf — 2 [9 ; namely that for low 
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temperatures T <C /i, /izj > /ig, implying that there is a 
phase with the thermodynamic properties of degenerate 
quark matter, but in which color is confined. The cur- 
rent result is if anything stronger than that of [5] because 
nq{^,) approaches the free quark result much more closely 
as a result of the smoother gauge fields. Such a phase 
is reminiscent of the confined, chirally-symmetric quarky- 
onic phase originally discussed in the context of large- iVc 
gauge theories [TB] . Unfortunately our use of Wilson lat- 
tice fermions, which have no chiral symmetry away from 
the limit k — Kc, precludes a discussion of whether chiral 
symmetry is restored for /i > fig at present. 
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Fig. 5. e,/e's"|j' vs. fj./m^ for Nf = 2,4. 



So far we have found a close similarity between Nf = 2 
and Nf = 4. This does not carry over to the quark energy 
density Sq, defined here by 

Nf 

^;(7o + l)e-^C/l_5^U)- (11) 

Fig. [5] plots eq/e'sB versus [i for Nf = 2,4. Note that a 
vacuum contribution evaluated at /i = must be sub- 
tracted from both interacting and free data; for N f — 4 
this correction e^a** — 0.3724(10). Even after this additive 
correction there is still a multiplicative renormalisation re- 
quired by a ^-independent factor known as a Karsch coef- 
ficient [T7] . Non-perturbative values for Karsch coefficients 
are still to be determined for QC2D, but the shapes of the 
curves are in principle correct up to discretisation errors. 
Fig. [5] shows a big difference at low values of ^ between 
N f — 2, where eq/e'^g has a peak considerably larger than 
that predicted by xPT [8p]. and Nf = 4, where the ra- 
tio actually starts negative and rises monotonically with 
[I. A negative value of Sq is not forbidden a ■priori, but 
the requirement for positivity of the total energy density 
certainly constrains the contribution Eg from the gluons 

(see below). For /i ?J /i_D eq'~'^ /s^sb ~ 2 becomes ap- 
proximately constant; e^^"''/^^^' approaches this value 



from below, and for /i/mjr ;C 1.5 the two models appear 
to coincide up to the unknown Karsch correction. 



.N=4 
N=2 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 
u/m 

Fig. 6. The unrenormalised quark contribution to the trace 
anomaly k,^^ Tr{l — (tp^p)) vs. /i/m,r for Nf — 2,4. 



A related quantity is the quark contribution to the 
trace of the stress-energy tensor {T^^)q, given by 

(T^^h = «7r ^ -i^NfN, - m). (12) 

With data from only one lattice spacing, we are currently 
unable to estimate the beta-function; Fig. [5] plots raw 
values of K~-'^Tr(l — {'4'4')) for Nf — 2,4, normalised to 
two quark flavors for ease of comparison, and including 
the necessary vacuum subtraction. Qualitatively they have 
very different behaviour for /i < and suggest that 
{Tfifi)q^~'^ and (T^^^)^^"^ differ even in sign in this regime. 
Note that xPT predicts T^^ > for /.to < /.t < ^3/io [S]- 
Since T^^ = e — 3p for isotropic matter, the negative sign 

of {Tfifi)q^~'^ is consistent with the negative value of eg for 
small /i reported in the previous paragraph. Once fJ. ^ hd, 
both models exhibit a strong upward trend, suggesting 
that quarks dominate T^^ in the deconfined phase. 

Finally we present results for local gluonic observables. 
In a non-Lorentz invariant system such as one with fJ- ^ 
it is helpful to define 



i<j 



3N, 



(13) 



We then consider in Fig.[7]the difference, proportional to 
the gluon energy density 



(14) 



where Z is another as yet undetermined Karsch coefficient 
(assumed unity in the figure), and in Fig. [5] the average 
plaquette related to the gluon component of the stress- 
energy tensor via 



(7)i^t)fl 



80 

-a/ X 3/?(n, 
oa 



□0 



(15) 
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Fig. 7. Eg/At"' vs. ^/m^ for Nf = 2,4. 



-0.05 



•-*Nf=4 
♦^♦Nj^2 



» 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 

u/m 

Fig. 8. r/ie unrenormalised gluon contribution to the trace 
anomaly 3/3(Ds + □«) fs. yL/m-n for Nf = 2, 4. 



For the gluon energy density (Fig. [7]) we have no an- 
alytic expectation from either free field theory or xPT. 
As in previous work, we plot the dimensionless quantity 
Eg/ fj,^; modulo the Karsch coefficient its value for ^ ^ /zq 
is approximately constant and similar in magnitude for 
Nf = 2,4 (at the smallest fi the errors are hard to con- 
trol), but for /i the ratio increases faster for Nf — 4, 
and has essentially doubled by fi/rriTr — 1-5. For (T^^)g 
(Fig. [8]) we note that again a vacuum subtraction must 

be applied. It is clear that there is a regime /J.o ~ M ^ M£> 
where (T^^)^^~^ is positive, before turning over to become 

negative at large fi; by contrast (T^^)^^"^ < through- 
out. It is not yet clear whether the sharp upward kink at 
/i/m^ = 1.3 is physical, or merely a discretisation artifact. 

The data of Figs. I6I8I taken together imply a major 
qualitative difference between Nf — 2 and Nf = 4: in 
the former case, T^^ > for all fi, which in particular 
is consistent with the existence of a non-relativistic de- 
scription with e ^ p for ^ ^ fio- For Nf = ihy contrast, 
T^n < for II < flu (we are unable to make a firm state- 
ment at this point about the deconfined phase due to the 



absence of the beta- functions), implying that the matter 
which forms at onset is already strongly-interacting and 
relativistic. 



4 Discussion 

In this paper we have presented the first exploratory re- 
sults for flavor-rich baryonic matter with Nf = 4 in QC2D. 
Since available resources have restricted us to a single 
lattice spacing, hopping parameter and diquark source 
strength, we have chosen to match the simulation to a pre- 
vious study of Nf = 2 on the same lattice size with the 
same value of m^n-a [S], implying the same onset chemical 
potential iJoO- measured in lattice units. Accordingly the 
comparison between the two theories is best performed us- 
ing observables measured in units of ttIt^. Since a^f^'^a? < 
cr'^f^^a'^, we expect the lattice gauge fields to be much 
smoother in the flavor-rich case, although that renders re- 
sults more susceptible to flnite- volume and thermal correc- 
tions due to the corresponding reduction in Ls\/a, Lty/a. 
One particular concern in the latter case is that the static 
potential V{R) used to calibrate the lattice could actu- 
ally be underestimated due to thermal screening. Another 
caveat which must eventually be addressed is control over 
the extrapolation to continuum a — > and zero source 
j — > hmits, currently being addressed for Nf = 2 [18]. 

Let us cite one example of a potential issue in the 
current study, the total energy density e. If we make the 
most optimistic assumption that the Karsch coefRcients 
are approximately unity, and that both j — > and a — > 
extrapolations will not change the results much, then the 
data of Figs. I5I7I suggest that for /i ~ /io with Nf = 4 
e//x^ = Eq/fi^ + Sg/fi^ ~ —0.7 + 0.13 < 0, which is clearly 
unphysical. 

We also note that at 11/111^^ = 1.76, the largest chemi- 
cal potential studied, the ratio Uq ^~^/n^^^ = 0.32, where 
^sut _ 2NfNc is the saturation density obtained when 
every lattice site is maximally occupied by fermionic de- 
grees of freedom; hence we might be concerned that in 
this high-density regime the results are susceptible to sat- 
uration artifacts. In this light, the superfluid condensate 
in Fig. m and gluon contribution to the stress-energy ten- 
sor in Fig. [5] are especially troubling, since both show a 
sharp kink at ^i/m^ = 1.3, for which we currently have no 
explanation. 

Bearing these caveats in mind, let us review the prin- 
cipal results. Fig. |3] for Uq, p for Nf = 2,4 expressed as 
fractions of the free-field results, demonstrates that both 
models display a qualitatively similar sequence of tran- 
sitions as fj, is increased. For Nf = 4 there is a range 
0.75 < n/niTr < 1.0 where riq/n}^^ « 1, implying the exis- 
tence of a well-defined Fermi surface. Further study will be 
needed to establish whether the difference in the height of 
this plateau results from the diminished impact of lattice 
artifacts in the Nf = 4 simulation, or a genuine physical 
difference between the two theories. Still, the enhanced 
flatness of the Nf — 4 plateau and the plausibility of the 
resulting physical picture increases our confidence in the 
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slightly ad hoc procedure we have used to compensate for 
discretisation artifacts. 

Another result is the confirmation of the major find- 
ing of [5], namely that hd > fJ-Q, or in other words, for 
both Nf = 2 and 4 there is a range of n where degener- 
ate quark matter remains color-confined. It will be very 
interesting to elucidate the nature of the deconfined phase 
further, since at first sight it bears little resemblance to 
thermal deconfinement. For instance. Fig. [7] suggests that 
the gluon energy density rises smoothly for /i ^ /xd, with 
Sg/fi'^ increasing by at most a factor of two over its value 
near onset. 

Finally, the most interesting and unexpected feature of 
the new simulation is the big difference between Nf = 2,4 
for the quark contributions to the energy density (Fig. [SJ 
and stress-energy tensor (Fig. For moderate fj, both 

eq^~'^ and (T^;^)^^"^ are negative; this result is self con- 
sistent and independent of the value of the relevant Karsch 
coefficient, though as discussed above could potentially al- 
ter as the limits a — 0, j — ?► are taken. The disparity 
between Nf = 2 and 4 however, and of each with the pre- 
dictions of xPT [5], is striking. It suggests very different 
physical descriptions in the low-density regime: as /i — > 
/Xo+ the Nf = 2 theory appears to be a non-relativistic 
BEC formed of weakly-interacting tightly-bound diquark 
bosons, consistent with xPT [14 and yielding T^^ > 0, 
whereas with Nf = 4 the matter appears to be relativistic 
and strongly- interacting for all fi > fio- 

This project was enabled with the assistance of IBM Deep 
Computing. S.K. was supported by the National Research Foun- 
dation of Korea grant funded by the Korea government (MEST) 
No. 2010-0022219. 



12. S.J. Hands and P. Kenny, [arXiv: 1104.05221 [hep-lat]] 

13. P. Kenny, Ph.D. thesis, Swansea University (2010). 

14. J.B. Kogut, M.A. Stephanov, D. Toublan, J.J.M. Ver- 
baarschot and A. Zhitnitsky, Nucl. Phys. B 582 (2000) 
477. 

15. S. Gupta, K. Huebner and O. Kaczmarek, Phys. Rev. D 
77 (2008) 034503. 

16. L. McLerran and R.D. Pisarski, Nucl. Phys. A 796 (2007) 
83. 

17. F. Karsch and I.O. Stamatescu, Phys. Lett. B 227, 153 
(1989). 

18. S. Hands and J.I. SkuUerud, in preparation. 



References 

1. S.J. Hands, I. Montvay, S.E. Morrison, M. Oevers, 
L. Scorzato and J.I. SkuUerud, Eur. Phys. J. C 17 (2000) 
285. 

2. J.B. Kogut, D.K. Sinclair, S.J. Hands and S.E. Morrison, 
Phys. Rev. D 64 (2001) 094505 

3. R. Aloisio, V. Azcoiti, G. Di Carlo, A. Galante and 
A. F. Grillo, Nucl. Phys. B 606 (2001) 322. 

4. S.J. Hands, I. Montvay, L. Scorzato and J. SkuUerud, Eur. 
Phys. J. C 22 (2001) 451. 

5. S.J. Hands, J.B. Kogut, M.P. Lombardo and S.E. Morri- 
son, Nucl. Phys. B 558 (1999) 327. 

6. S. Hands, P. Sitch and J.I. SkuUerud, Phys. Lett. B 662, 
405 (2008). 

7. S. Muroya, A. Nakamura and C. Nonaka, Phys. Lett. B 
551 (2003) 305. 

8. S. Hands, S. Kim and J.I. SkuUerud, Eur. Phys. J. C 48 
(2006) 193. 

9. S. Hands, S. Kim and J.I. SkuUerud, Phys. Rev. D 81 
(2010) 091502(R). 

10. M.P. Lombardo, M.L. Paciello, S. Petrarca and B. Tagli- 
enti, Eur. Phys. J. C 58, 69 (2008). 

11. B. AUes, M. D'EUa and M.P. Lombardo, Nucl. Phys. B 
752, 124 (2006). 



